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Abstract 

Acoustic receptivity of a Blasius boundary layer in the presence of distributed surface irregular- 
ities is investigated analytically. It is shown that, out of the entire spatial spectrum of the surface 
irregularities, only a small band of Fourier components can lead to an efficient conversion of the 
acoustic input at any given frequency to an unstable eigenmode of the boundary layer flow. The 
location, and width, of this most receptive band of wavenumbers corresponds to a relative detuning 
of 0(R~^ S ) witli respect to the lower-neutral instability wavenumber at the frequency under con- 
sideration, flu, being the Reynolds number based on a typical boundary-layer thickness at the lower 
branch of the neutral stability curve. Surface imperfections in the form of discrete mode waviness 
in this range of wavenumbers lead to initial instability amplitudes which are 0(7? 3 ^ 8 ) larger than 
those caused by a single, isolated roughness element. In contrast, irregularities with a continuous 
spatial spectrum produce much smaller instability amplitudes, even compared to the isolated case, 
since the increase due to the resonant nature of the response is more than compensated for by the 
asymptotically small band-width of the receptivity process. Analytical expressions for the maximum 
possible instability amplitudes, as well as their expectation for an ensemble of statistically irregular 
surfaces with random phase distributions, are also presented. 



1 Introduction 


The purpose behind this paper is to present some theoretical results concerning the acoustic receptivity 
of a boundary layer flow due to distributed surface imperfections, which can be either deterministic 
or random in terms of their origin. Such imperfections may arise due to a variety of causes, ranging 
from manufacturing defects or structural joints to operational factors, such as paint erosion, insect 
debris and ice accretion, etc. [1). Of the many ways they have been known to afTecl the transition to 
turbulence, inducement of receptivity is a major one when the maximum height of these irregularities 
is small. 

The generation of Tollmien-Schlichting (henceforth T-S) instabilities in low-speed boundary layers 
via the interaction of free-stream sound with a local surface distortion was first explained by Goldstein 
[2] and Ruban [3] in 1985. They showed that receptivity in this case can be attributed to the fact 
that the unsteady scattered field produced by this interaction inherits its temporal scale from the 
free-stream disturbance but its spatial scales from the shorter surface irregularity, thereby acquiring a 
Fourier spectrum that overlaps with the T-S wave. Following this fundamental breakthrough, a variety 
of other localized receptivity problems were studied by other investigators, and examples of these can 
be found in the proceedings [4] and [5]. 

Receptivity due to distributed waviness of the surface was studied by Crouch [G] [7]; however, the 
conclusions of these two references concerning the magnitude of receptivity are contradictory to each 
other. Furthermore, the prediction technique used therein appears to be more suitable for distributions 
of surface nonuniformities which are homogeneous in the flow direction and only involve a finite number 
of (discrete) Fourier modes. On the other hand, Choudhari and Street! [8] had shown how generation 
of instabilities via distributed regions of surface inhomogeneities could, in general, be predicted by a 
relatively simple extension of the localized receptivity results. Ibis extension is based on treating the 
receptivity in a region of large streamwise extent as a sum of the contributions from its constituent sub- 
domains of infinitesimal size, in each of which localized receptivity analysis becomes applicable. The 
same idea was used originally by Tam [9] in an attempt to predict the direct excitation of instability 
waves in the absence of any non-uniformities on the surface. The advantages of this technique include 
its applicability to arbitrary distributions of surface inhomogeneities, and its simplicity in leading to a 
first order ordinary differential equation for the amplitude of the instability wave at any given frequency. 
The latter, as demonstrated in this paper, facilitates a further analytical treatment of the distributed 
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receptivity problem, thereby yielding closed form expressions for the generated instability motion and 
also providing a clearer interpretation of the physics involved. 

The receptivity to a single frequency acoustic disturbance in the presence of surface imperfections 
of known shape and distribution is considered first in Section 2.2 below, following a brief discussion of 
the general procedure for solving receptivity problems involving weak surface inhomogeneities; a more 
detailed discussion of the background analysis can be found in Refs. [8] and [10]. Although the classical 
Orr-Sommerfeld framework has been utilized as the basis throughout, a similar treatment could also 
have been applied to the problem using a purely asymptotic theory. The results obtained in Section 
2.2 for the deterministic problem are used as the building block in Section 2.3 for analyzing receptivity 
due to a random roughness distribution with a specified power spectral density. Numerical verification 
of the theoretical results is also considered. 

2 Receptivity Analysis for Distributed Surface Inhomogeneities 

2.1 Background and Problem Formulation 

Consider the two-dimensional flow over a slightly rough surface that is nominally flat and aligned with 
an incompressible free-stream with mean speed plus time-harmonic fluctuations having frequency 
w*, and amplitude u* c , such that ej 3 = u* c /U ^ , << 1 (Fig. 1). The maximum perturbation in the 
surface height x *) with respect to the mean position y* — 0 is also assumed to be sufficiently small 
compared to the mean boundary-layer thickness (c w = (/i* ; /<5*) ma:r << 1) so that the (low everywhere 
can be expanded as a regular perturbation series in terms of the amplitude parameters e w and ej s . As 
described below, each term in this perturbation series turns out to involve a different combination of 
spatial and temporal scales, and the objective from the standpoint of receptivity is then to determine 
the first term which contains the desired combination of scales, i.e., one that overlaps with the local 
instability wave. Of course, in order to determine this term with a given order of accuracy in terms of 
Reynolds number effects, one needs to know all the previous terms also to the same order of accuracy. 

It should be obvious that the zeroth order term corresponds to the Dlasius streamfunction for 
the mean boundary-layer flow over a flat surface, while the 0(cj s ) perturbation, corresponding to 
the unsteady signature of the acoustic frec-strcam fluctuation within this boundary layer, is given by 
the generalized Stokes-wave solution obtained in [11]. The 0(c w ) perturbation, corresponding to the 
steady but short-scale disturbance due to variations in the surface geometry, satisfies the parallel flow 
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equations up to 0(R~ 3 / 4 log R) on a local basis, although the large-scale modulation due to the weak 
growth of the Dlasius boundary layer also comes into play, depending on the streamwise extent of the 
region of imperfections. Note that R denotes the Reynolds number based on the boundary-layer scale 
and is defined as 72 = \/Re x * y Re x ♦ being the Reynolds number based on the distance from the leading 
edge. 

The first term in the amplitude expansion which possesses both unsteadiness and short spatial scales, 
and hence is the object of the receptivity calculation, corresponds to the 0(e w cj 9 ) term arising from 
a quadratic coupling between the unsteady free-stream and steady surface disturbances. Physically, 
it represents the unsteady field produced by the interaction of the 0(ej 3 ) Stokes shear wave with the 
0{c w ) mean flow perturbation, and by its direct scattering due to the 0{c w ) perturbation in the surface 
height. Within the “geometrical optics 1 ’ approximation, this scattered field is governed by the unsteady 
linear disturbance equations for a locally parallel mean flow, along with inhomogeneous terms in the 
differential equations themselves and in the boundary conditions as well, in view of the two different 
mechannisms for scattering as mentioned above. Again, non-parallel effects will need to be considered 
if propagation distances become comparable to the body-length scale. However, subsequent analysis 
will show that most of the receptivity is actually concentrated in a region of much smaller streamwise 
dimension and hence, a prediction based on the parallel flow equations is sufficient, at least to the 
leading order of approximation in terms of the Reynolds number. 

Doth the mean flow perturbation and the unsteady scattered field, thus, satisfy the quasi- parallel 
disturbance equations, and therefore can be solved on a local basis by reducing these partial differential 
equations to ordinary ones through the use of a Fourier transform in the streamwise direction. The 
equations in the transform space can be solved either in closed form, via a systematic but multi- 
layer expansion in terms of the Reynolds number [2-5], or numerically by using the Orr-Sommerfeld 
equation as a non-asymptotic but composite approximation [4-8] . In either case, the locally generated 
instability wave can be isolated as the residue contribution corresponding to a pole singularity of the 
inverse Fourier integral for the unsteady scattered field. This residue contribution is linear in terms of 
the local perturbation to the mean surface-height, with a coefficient function, A u , that characterizes 
the efficiency of the local receptivity process. The T-S fluctuations produced locally, i.e., in each 
infinitesimal subregion, propagate independently of the instability motion generated elsewhere, with 
their amplit ude and phase variation being determined by the T-S wavenumber as a function of R. The 
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total instability amplitude at any location is then given by the integral over the contributions from all 
sources upstream of this locatioiJ 8,10 ^ 


/ _ R 

- E U {Y ; R) e’t©r5(/i)-^J f h w (R a ) c"' 9 ^ ^ II (R - R s ) dll 3 , (2.1a) 

7T J 


where u* TS denotes the streamwise velocity fluctuation associated with the T-S motion, E u is the T-S 
eigenfunction for this quantity, normalized to have a maximum magnitude of unity across the boundary 
layer, and Qrs is the spatial phase of the T-S wave 


Qrs(R) = 2 jc* TS (R) dll . 


( 2 . 16 ) 


The quantity h w (x *) is the local perturbation in the surface height, nondimensionalized by the slowly 
varying length scale x*Re x . , and A u is the efficiency function obtained from localized receptivity 
analysis. One should note that the streamwise variations of the eigenfunction E u , wavenumber qts> 
and efficiency factor A u are much slower than those of the instability phase Qts and the surface-height 
function h w ; thus we have | << 1 etc. In Fig. 2, we have shown the Reynolds number 

dependence of the magnitude of A u for selected values of the nondimensional frequency parameter 


/o = 10 6 x u)*v* /U^. 

By differentiating with respect to It and invoking the quasi-parallel approximation, Eq. (2.1a) can 
also be converted to a wave amplitude equation similar to that obtained by Tam [9], 


d u TS 

(lit 


~ 1 


. </0 


TS 


(lit 



(2.2) 


where, for simplicity, we have omitted the Independence of the disturbance motion, and from now 
on, u? rs will denote simply the maximum of the streamwise velocity perturbation at each streamwise 
location. The homogeneous solution to this amplitude equation corresponds to the T-S wave, which is 
excited every time the local spectrum of the geometry function h w overlaps with the T-S wavenumber. 
As seen from Eq. (2.2), when the receptivity occurs continuously over a large number of instability 
wavelengths, the change in the wave amplitude at any station is a combined outcome of t lie local 
amplification of the instability waves generated upstream of the present location and the external input 
due to local receptivity. As discussed in Refs. [9] and [10], the external input dominates the initial 
development of the instability amplitude; however, after u^ s becomes sufficiently large in magnitude, 
the amplitude evolution curve asymptotes to that of a pure T-S wave eigensolution. Note also that, 
as a result of the distributed receptivity process, the location of stationary (i.e., maximum) amplitude 
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does not correspond to the theoretical upper-branch location R u , b .. However, the shift is asymptotically 
small, and one may, therefore, measure the generated instability motion in terms of its amplitude at the 
upper-branch location itself. The intrinsic receptivity can be gauged by dividing the nondimensional 
amplitude u* rs {R = Ru.b.)/Kc by the amplification ratio between the two neutral locations, which 
leads to an “effective coupling coefficient”, C, indicating the effective instability amplitude at the 
beginning of the linear amplification stage. A coupling coefficient^ basically relates the output of the 
receptivity process, i.e., the amplitude of the generated instability wave, to its input, i.e., the free- 
stream disturbance amplitude, and the measure C is the local coupling coefficient for an equivalent, 
but ficticious, localized mechanism that lias all its receptivity lumped together at the lower-branch 
station R = Ri. b .. Approximate analytical expressions for the effective coupling coefficient for different 
types of geometries arc derived in Sections 2.2 and 2.3 below. 

2.2 Receptivity due to a single-mode surface waviness 

Since any spatially homogeneous distribution of roughness elements can be examined in terms of its 
individual Fourier components, we first consider geometries of the form 

h w (R) = J -^e iQ ^ n \ where Q w = I* 2a w {R)dR + fa, a w (R) = , (2.3) 

corresponding to a surface with single-mode waviness of constant nondimensional amplitude Rh 0 - Note 
that Rh 0 is the Reynolds number based on the free-stream speed U and the dimensional amplitude 
of the surface-height perturbation Iiq. Similarly, R,\^ is the Reynolds number based on the constant 
dimensional wavelength A£,(= 27 r/o^) of the surface undulations, while a w denotes the dimensionless 
wavenumber based on the slowly-varying, local length scale, x*Re x . . The integral solution (2.1a) 
now becomes 

n 

u 'TS _ [2. Ji h e ie TS (rt)+i<t>o-ivt f e i[0u.(R.)- 0 rs(n»)] ( l]{ s (2.4) 

u* ac V 7T 0 J Rs 

Since the efficiency function A u varies slowly with the integration variable R s , the integrand in (2.4) 
is nearly oscillatory if the wavenumber of the surface undulations differs from the local instability 
wavenumber ay s by an 0(1) magnitude. Then contributions from the neighbouring surface locations 
corresponding to the out-of-phase elements during each cycle of oscillation nearly cancel each other, and 
the integral becomes end-point dominated. This implies that the disturbance amplitudes will retain 
the the same order of magnitude throughout the region of waviness. 
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2.2.1 Waviness synchronized with the neutral instability motion 

A sustained generation of instabilities will occur only if n w « c*rs> corresponding to a nearly syn- 
chronous variation of the forcing phase 0 U; and the eigen mode phase Ojs-> which leads to a mutual 
reinforcement of the contributions to the integral in (2.4) from adjacent surface locations. Of course, 
since the instability wavenumber varies slowly with distance along the surface, such a resonance is 
necessarily localized in space, and the extent of this localization determines how large the distributed 
receptivity is, when compared to the receptivity induced by a single, isolated roughness element. Since 
the instability wavenumber is complex, in general, the condition of perfect local resonance requires that 

«; = a* rs {R,. b .) or a* w = a’ rs (R u . b .) , (2.5a, h) 

where the asterisk denotes dimensional quantities, and subscripts Lb . and u.fc. represent lower and upper 
branches of the neutral stability curve, respectively. When (2.5a.) or (2.5b) is satisfied, the exponent 
in the integrand of (2.4) has a saddle point at the respective neutral location, with the real R axis 
corresponding to a path of descent in case of (2.5a), and a path of ascent when (2.5b) is true. The 
latter case is, of course, of little practical interest, and will not be discussed here. 

The presence of a saddle point at R — Ri.b, implies that the importance of the local contribution to 
the effective coupling coefficient diminishes rapidly with distance away from the lower branch location. 
Applying the steepest descent method yields the following approximation for the effective coupling 
coefficient 

Cp.s. = -Jr 2 - t/r|-A u (7?,.6.)e^ 0 , D a = a' w - a' TS , (2.6) 

Ri.b. V li) <* 

where the suffix p.s. indicates that the condition (2.5a.) for perfect synchronization is satisfied, and 
the primes denote derivatives with respect to R , evaluated at the neutral location R = 72/6.. The 
dispersion factor D a , which is simply in dimensional terms, controls the rate at 

which the phase of the natural response becomes uncorrelated with the surface geometry as one moves 
away from the resonance location. Asymptotically, the magnitude of this factor is of 0(/Z| b ^ 4 ), and 
therefore the length of the resonance region corresponds to AT? = 0(7 ?j 5 { 8 ), or in physical terms, 
Ax* = 0(x* b Tij”^ 8 ), which is of the same order as the geometric mean of the distance from the 
leading edge and the instability wavelength at the lower branch location. This shows that that the 
synchronized surface waviness will lead to instability amplitudes which are larger than those produced 
by an isolated surface protuberance by a factor of This simple analytical result confirms, 
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as well as explains, the more recent computations of Crouch (Ref. 7) where he found the distributed 
receptivity to be “two orders of magnitude” larger than that due to an isolated roughness element, as 
against his earlier findings (Ref. 6) that the receptivity is comparable in both cases. Finally, since 
in practice, the imaginary part of D a is larger than its real part at all frequencies of interest, the 
dispersion effect which leads to the detuning of the forcing 0^(1?) with respect to Ots(R) caused 
more by the variation in the instability growth-rate, than by a variation in the instability phase-speed 


near R = Ri,b., 

It is equally instructive to consider the single mode waviness problem using Tain’s wave amplitude 
equation (2.2), which is physically analogous to the second order equation obtained for a forced oscillator 
with slowly- varying parameters, except that the oscillator motion usually involves a zero or positive 
damping. Both linear and nonlinear problems of this type have received, and continue to receive, 
ample attention from applied mathematicians over the years as evidenced from the various citations in 
a recent review by Kevorkian [12], which appea.r to have begun with Itef. [13] in 1971. Hence, only a 
brief sketch of the solution to (2.2) is given here. A reader interested only in receptivity may also find 
some supplemental discussion, especially related to the appropriate initial condition for (2.2), as well 
as the numerical results for the case of wall-suction induced receptivity, in Ref. [10]. 

Basically, the total disturbance amplitude consists of a superposition of the particular solution, 


U p = U ac 


1 Rh 


A«(*i.6.) 


'MA'o.ft.y e z©^(/?) 

i\/27r Rib. a w (Ri,b,) - ors{Kl.b.) 


and a constant multiple of the homogeneous solution, 


(2.7a) 


I 


i 


t* = n* p*®Ts{R) 
A ac*" 


(2-7 b) 


If (\ w /ajs — I = 0(1) everywhere, then the two solutions (2.7a) and (2.7b) remain decoupled, i.e., 
linearly independent, everywhere, and imposing an initial condition corresponding to zero initial in- 
stability amplitude forces the coefficient of the homogeneous solution to be zero identically. On the 
other hand, when the resonance condition (2.5a) is satisfied, the particular solution becomes singular 
at R = Rib., reflecting the fact that the particular and homogeneous solutions are linearly dependent 
in the vicinity of this location, and, therefore, the decomposition of the total solution into (2.7a) and 
(2.7b) is invalid. 

_ ^ n to 

A separate, inner expansion becomes necessary for (/£ — Rib.)/ ftl.b. ~ 0(R lb ), wherein both a w (R) 
and r*7\s(7£) can be approximated by their respective linear Taylor-serics approximations centered on 
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the resonance location, R — Ri,b.- Expressing the coupling coefficient Uj S /u* ac in the inner region as the 
product of a rapidly-varying component e 2tQl b ^~ Rl b \ which incorporates the phase synchronization 
with the local instability motion, with an amplitude function A(R) that varies on the inner length-scale, 
we obtain the amplitude evolution equation 


dA 

dR 


= 2 ia' TS (R - Ri.b.)A + a’ w A u (R, b .) e i *°+ i <( R - R < »y, (2.8) 


where, for brevity, we have avoided introducing a specific coordinate just for the inner region. The 
solution to (2.8), 


A(R) 


A,. b e ia Ts( n - R *■*•>’ + K(ll,. b .) c iQ Ts( R - R ‘ »y [ R j*o+iD a (R,-n l b ? dIis J (2 . 9) 

V 7T Ri k J Hi h 


Rl.b . JRlb. 

matches the particular solution (2.7a) sufficiently far upstream, provided 


yj ^*0 A u ( Rl.b . ) ^i(j, 0 

,b ■ ~ Rl.b. y/2 iD~ a 


( 2 . 10 ) 


As R — Ri. b . —* oo on the inner scale, Eq. (2.10) yields 


u* TS /u ac -► 2 Ai. b . e 2ia Ts(iii. b Mn-n t b .)+io,' TS (n-n,. b ) 3 
_ A|i(/?;.fc ) e i<t>o+2ia'Ts(Ri b .)(n-Ri.b ) 2 (2 11) 

Rl.b. is/2^ D a (R- Rl.b.) 

implying that the disturbance motion downstream of the inner region is dominated by the homogeneous 
(i.e., the instability wave) solution, which is larger in amplitude than the particular solution (2.7a) by 
a factor of at least (R - Ri. b .)y/D^, = 0(7if{ 8 ); see Fig. 3. The effective coupling coefficient in the 
perfectly synchronous case is then given by 

C P .M ) = 2/ln,. = A u (R,.b.) e*\ (2.12) 

which is identical to the steepest-descent result (2.6). The reader may note that, since the coupling 
coefficient was determined entirely by the solution in the shorter inner region, nonparallel effects asso- 
ciated with the streamwise derivatives of the base flow, and the vertical velocity associated with the 
same, can be neglected to the leading order, at least. 


2.2.2 Near resonant geometries: the effect of detuning 

Other types of near- resonant geometries, corresponding to arbitrary modulation of the waviness am- 
plitude on the length scale of the resonance region, can also be handled in a manner similar to that 
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in Section 2.2.1 above. Thus, for geometries of the form h w = A w (R)e 2 ' ars ^ R,b ^ R n,li \ with A W (R) 
varying on the inner scale, one obtains 


1(7 i) = /l,.i.c' v ’rs (n - ,?lk )2 + t/? A u{Rl.k.) e ia r*< n /? "’ )2 [ A w (R,)e 

V 7T it/, b. JIh 6, 


m' *(/?,-/*! 


<//f s . (2.13) 


For the case of spatially periodic geometries with a wavenumber detuning of a wi with respect to 
otTs(Rt.b.), i-e., 

A W (R) = e 2ioru 'i {/?-ni6 ) +, ^o, (2.14a) 


this leads to 

C(a m ) = € *Ua C p . s ., 


(2.14h) 


indicating that the effective coupling coefficient due to an individual Fourier component drops off in 
a gaussian manner away from the resonant wavenumber Q/.t.. Ily observing that only the real part of 
iD a contributes to the magnitude of the exponential factor, it is clear that this drop-off is purely due 
to the streamwise variation in the instability growth rale near R = Rib.- Basically, if one views the 
phase-synchronization process in terms of the real parts of the respective wavenumbers [a w and Q75), 
then, for the detuned cases, the center of the resonance region is somewhat upstream or downstream 
of the lower branch location, which implies a reduced amplification ratio between the resonance region 
and the upper branch, and in turn, a reduced magnitude of the effective coupling coefficient. Only 
a detuning of up to = 0{\/TJ^), corresponding to a w JaTs{Ri.b.) = 0(R t ^ S ), can produce any 
significant receptivity, due to the scaling of A/Z/TZi.t. = 0(/? ( 6 3 ^ 8 ) for the resonance region. 

Figure 4 shows the magnitude of the normalized effective coupling coefficient, C(a Wi )/ Rh 0 , plotted 
as a function of the wavenumber detuning parameter ct Wl /<*Ts(Ri.b.) expressed as a percentage, for two 
different values of the frequency parameter / = u*u*/U^. The frequency of / = 25 x 10 G has a loga- 
rithmic amplification factor of about 9, and, hence, lies in the most critical range of frequencies for the 
transition process. On the other hand, the frequency of / = 55 x 10 -6 , although not very important for 
transition, may be relevant to laboratory experiments. The figure shows that the analytical predictions 
at both frequencies are in very good agreement with the results obtained from a numerical solution of 
the wave amplitude equation (2.2). The only significant difference between the two results appears to 
be that, at large positive values of the detuning parameter (> 8%), the numerical coupling coefficient 
decreases more slowly than the gaussian drop-off predicted in the analytical solution (2.1 lb). 


10 


1 Mu mill miiiiiiiiii 



2.3 Geometries with a continuous spectrum of a deterministic or random nature 

Now, consider the class of surface imperfections which have a continuous spatial spectrum, and, there- 
fore, cannot be classified as discrete-mode waviness. Assuming the spectrum to be homogeneous in the 
streamwise direction, one can integrate over all contributions of the form (2.14b) in order to arrive at 
the total coupling coefficient for any particular geometry. In addition to its dependence on the power 
spectral density R\ Q (a Wl ) of the surface-irregularity distribution, the coupling coefficient also depends 
upon the relative phase distribution, <f> { 0 (a Wl ), between the different wavenumber components of the 
surface undulations. However, an obvious conclusion that holds for all geometries with a continuous 
spatial spectrum is that the integrated coupling coefficient is 0(R t ^f S ) times smaller than the discrete 
mode solutions (2.12) and (2.14a), since the bandwidth in the space which can contribute signifi- 
cantly to this integral is only of 0(RJ’^ 8 ). This makes the overall receptivity even smaller than that 
induced by an isolated roughness element with a comparable perturbation in height. Note that the 
receptivity is stronger for the discrete inode waviness of the form (2.3), because the latter represents 
a peak in the spatial spectrum of the surface irregularity near the resonant wavenumber, and, thus, 
possescs more energy in the relevant band of wavenumbers. 

Lacking a precise knowledge about the surface irregularities, it is not an easy matter to characterize 
the phase distribution of the spectral representation of any stationary process. Only a few limited results 
exist for finite-dimensional distributions [14]; these prove that the phase at any wavenumber is always 
uniformly distributed on [ — 7r, 7 t] and that the joint distributions of the phase variables are independent 
of the spectral amplitudes. Assuming that this result also holds in the limiting case of a continuous 
spectrum, and that the corresponding power spectral density R\ o does not vary significantly over the 
narrow bandwidth of the spectrum, one finds that the expected value < C 2 > for an ensemble of 
rough surfaces which are statistically similar is given by 

< C 2 > 1/2 = sin(arg(D a )) A u (#/.6.)| , (2.15) 

Ri.b. 

where arg denotes the argument of a complex valued quantity. In any given realization, the effective 
coupling coefficient will vary from a minimum value of zero to a maximum of 2 3 / 1 times the average 
value in (2.15). This upper limit is associated with that particular realization for which contributions 
from all wavenumbers are in phase, and, therefore, the respective magnitudes can be added to each 
other in a linear fashion. 
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3 Concluding Remarks 

The theory presented in this paper elucidates the nature of the distributed receptivity process, and 
provides a closed form approximation for the instability-wave amplitude as a function of the amplitude 
and frequency of the free-stream disturbance. Although analyzed in the specific context of receptivity 
due to distributed surface irregularities, the results obtained can be easily generalized to receptivity 
due to other forms of surface disturbances, such as wall suction, wall admittance/compliance, and wall 
temperature, etc. An especially important application would be to predict the generation of stationary 
cross-flow vortices due to distributed surface roughness/suction by generalizing the localized receptivity 
results obtained earlier [15]. 

The theory also demonstrates the increased effectiveness of near-resonant forcing, which leads to 
continual receptivity over long distances. At low speeds, such resonant forcing can only be produced 
via the interaction of the free-stream disturbances with surface inhomogeneities of different types, such 
as surface roughness, short-scale variations in suctions, etc. However, in supersonic flows, phase speeds 
of the free-stream and boundary-layer disturbances are of the same order, and a near-resonant forcing 
could possibly arise directly, without any need for “tuning” via scattering at the boundary surfaces. If 
such an unsteady forcing does indeed occur, the ideas presented above might turn out to be applicable 
in predicting the boundary-layer response in these problems as well. Our current effort is focussed 
along this direction. 
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Magnitude of the efficiency factor \ from Eq. (2.1a) 
as a function of the vail inhomogeneity location R 

with / 0 = 10 6 x wV/iC as a parameter. 

The lower- and upper-branch neutral locations at each frequency 
are indicated by open and filled circles, respectively. 




Fig. 3 A sketch of the development of the wave amplitude across the 
resonance region. Note that the decomposition of the solution 
between the particular and homogeneous components inside the 
resonance region is purely imaginary, and has been shown for 
the purpose of intuitive understanding only. 


17 



\C\/R 



Fig. 4 Receptivity induced by discrete-mode surface vaviness: 
the normalized effective coupling coefficient, \C\/R^, 

as a function of the vavenumber detuning parameter Aa=a Wi /a rs (/? Z6 )xl00 
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